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Abstract 

It is beneficial to use a distributed parameter model for large space 
structures, because the approach minimizes the number of model 
parameters. Holzer's transfer matrix method provides a useful means to 
simplify and standardize the procedure for solving the system of partial 
differential equations. Any large space structures can be broken down into 
sub-structures with simple elastic and dynamical properties. For each single 
element, such as beam, tether, or rigid body, we can derive the 
corresponding transfer matrix. Combining these elements' matrices enables 
the solution of the global system equations. The characteristic equation can 
then be formed by satisfying the appropriate boundary conditions. Then 
natural frequencies and mode shapes can be determined by searching the 
roots of the characteristic equation at frequencies within the range of 
interest. This paper applies this methodology, and the maximum 
likelihood estimation method, to refine the modal characteristics of the 
NASA Mini-Mast Truss by successively matching the theoretical response to 
the test data of the truss. The method is being applied to more complex 
configurations. 


1. Introduction 

Control of flexible spacecraft is best analyzed by representing, in a 
single set of equations, all of the structural modes and the control system 
dynamics. Distributed parameter models enable this approach based upon 
the classical partial differential equation theories. In the early 1960's, the 
distributed parameter approach was developed simultaneously with the 
lumped parameter approach [1-3]. With the event of high-speed and large- 
memory computers the finite element method has been developed much 
more extensively than the distributed parameter approach. However, the 
advantages of the distributed parameter approach to control synthesis, 
parameter estimation and integrated design have been largely neglected by 
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the technical community. It is the purpose of this paper to show the 
advantage of employing the distributed parameter approach to large space 

structures. . . 

Holzer's transfer matrix method provides a useful means to simplify 
and standardize the procedure for solving the system of partial differential 
equations. Also, the transfer matrix method enables the construction of a 
relatively simple mathematical model for complicated structures. It is of 
great practical value to take advantage of catalogs of the most important 
transfer matrices readily available. A complex structure, then, only requires 
combining individual matrices to represent the structure which consists of 
connected elements. The similarity of this concept to that of the transfer 
function is particularly useful to control analysis. 

The transfer matrix method itself has been a matured method [4], and 
jfg power has been shown in several technical areas [5,6]. But, little use of 
the transfer matrix method has been made for the distributed modeling, 
parameter estimation and control of large flexible space structures. 

In this paper we applied the transfer matrix method, accompanied by 
the maximum likelihood estimation technique, to estimate the lateral 
bending characteristics of the NASA Mini-Mast truss (Fig.l. l)[7] by matching 
the theoretical transient response to test data. The Mini-Mast truss is a 
ground testbed for the Control-Structure Interaction (CSI) program. The 
total height of the truss is 20.16 meters, containing 18 deployable bays. 
Two instrumentation platforms have been installed at Bay 10 and Bay 18. 
Mini-Mast has 162 major structural elements. Finite element models of the 
truss involve thousands of elements. The distributed parameter model of 
the Mini-Mast truss used in this paper consists of two flexible beam elements 
and two rigid bodies. 

The shear deformation of the truss requires a Timoshenko beam 
model in order to match the frequencies at higher mode numbers. It is also 
necessary to extend the simple model by adding the so-called "appendage 
model" to account for the effects of dynamics of diagonal struts and 
associated hinge bodies. The method of this paper is shown to be applicable 
to more complex configurations. 


2. Derivation of Transfer Matrix 

Holzer's transfer matrix method [4] provides a useful means to simplify 
and standardize the procedure for solving the partial differential equations. 
Any large space structures can be broken down into sub-structures with 
simple elastic and dynamical properties. For each single element, such as 
beam in bending, rigid body, we can derive the corresponding field matrix 
and point matrix. Combining these elements' matrices in a required 
manner, one can calculate the responses, i.e. the solution to the global 
system equation, by proceeding from one point of the system to the other. 

This paper concentrates on the estimation of the lateral bending 
frequency of the NASA Mini-Mast Truss. The truss is modeled as two 
successive beam elements with two rigid bodies at the Bay 10, and Bay 18 (tip 
of the truss), respectively. To derive the transfer matrix we consider a 
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^ 1.1 NASA Mini-Mast Truss 

“original photo not available’ 


79 



cascaded beam-body system (Fig.2.1). The jth section consists o f a flexible 
beam whose elastic and dynamical properties will be described by a fiel 
matrix, and a rigid body whose dynamical property is presented by a point 
matrix. When necessary, we will designate any corresponding quantities to 
the left and right of a rigid body by superscripts L and R. 



Fig. 2.1 A Cascaded Beam-Body System 


The lateral bending of the beam is represented by the Bernoulli-Euler 
beam equation, 

-x 4 -\2 

+ = o (2.1) 

dz 4 a 2 at 2 

where, y(z,t) is the lateral displacement, and a 2 =k/m where k=EI bending 
stiffness and m=pA mass per length of the beam. By separation of variables 
y(z,t)=Y(z)T(t), we have two ordinary differential equations in Y(z) and T(t), 

Y""(z) - p 4 Y(z) = 0 ( 2 - 2 ) 

and 

T(t) + to 2 T(t) = 0 ( 2 - 3 ) 

where, p4=co 2 /a 2 . The solution to the Eq.(2.2) has the form, 

Y(z) = AsinPz + BcosPz + CsinhPz + DcoshPz (2.4) 

At the left end of the beam (z=0), the displacement Y(0), slope Y'(0), shear 
Q(0), and bending moment M(0) will be 
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*fl - Bj + Dj 

Yf l = Ajft + Cjpj 

Qf-l - -kjY/"(0) - Ajkjpj ■ CjkjP? 

M ?-l - -kjYj"(0) - BjkjpJ - Djkjpj 


or written in matrix form as, 
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Similarly, at the right end of the beam (z=Lj), the corresponding quantities 
are, 


Y j L = AjsinpjLj + BjCospjLj + CjSinhpjLj + DjcoshpjLj 

y ’j L = AjPjCOSpjLj - BjPjsinPjLj + CjPjCoshPjLj + DjPjsinhPjLj 

Q^=kYj , "(L)=-Ajkjp 3 cosp j Lj + B J k J p 3 sinp J L J + CjkjpJcoshpjLj + Djkjp-sinhPjLj 
M[=kYj"(L)=-A j k J p J 2 sinp J L j - BjkjpJcospjLj + CjkjpfsjnhpjLj + DjkjPjCOshpjLj 

or written in matrix form as, 

sinPL cospL sinhpL coshpL 

PcosPL -PsinPL PcoshpL PsinhpL 

-kp 3 cosPL kp 3 sinpL kp 3 coshpL kp 3 sinhpL 
-kp 2 sinpL -kp 2 cospL kp 2 sinhpL kp 2 coshpL 


b A 1 


( 2 . 6 ) 



Substituting Eq.(2.5) into Eq.(2.6) we obtain 
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YIL |Y\k 

H - [nil j Y ' 

Q ’ Q 

Ml i * M J j.i 


where, the field matrix of the jth beam element. 


[FM] j = 


l<cosBL+coshpL) -ksinpL+sinhpL) -±-<sinPL-sinhpL) -L(cospL-coshpL) 

r 2p 2kp 3 2kp 2 

-£<sinpL-sinhpL) l<cospL+coshpL) -i-(cosPL-coshpL) - -L<sinpL+sinhpL) 

2 2 2kp 2 •“P 

ikp\sinpL+sinhpL) - !kp 2 (cospL-coshpL) - l<cosPL+coshpL) &{sinpL-sinhpL) 

2 2 2 l 

-ikp 2 (cospL-coshpL) - Ikp(smpL-sinhpL) - -UsinpL+stnhpL) - l<cospL+coshpL) 

2 2 2p l 



Fig. 2.2 Free Diagram of a Rigid Body 

Next, let us consider the jth body (Fig.2.2). The translational and 
rotational motions of the body can be described by the following equations, 


mj Vcm = Of - Qj' ( 2 - 8 ) 

Ij y' L = Mf - Mj' - Qfrf - Qfif ( 2 . 9 ) 

For homogeneous motions, Eqs.(2.8) and (2.9) can be written as, 

Qf = Qf-m j O) 2 Y cm (2.10) 

M r = Mj L + Qfrf + Qfrf - IjO> 2 Y'f (2. 1 1) 


But, the displacement of the center of mass, Y cm . is related to Yj R and Yj L by, 
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From Eq.(2. 13), we see that 


( 2 . 12 ) 

(2.13) 


yR _ y r R yiR 
l ) ~ icm Ij I j 

Y L = Y + X y’L 

1 j 1 cm ^ I j I j 


Y — _ t 4 " Y ,l 

icm ’ Ij IJ I j 

To keep the compatibility of the deflection, it must be that 

ytR yiL 

J “ 1 j 

Substitution of Eqs.(2.14) and (2.15) into Eq.(2.12) gives 

Y R — V L T* V'L 

V, L P J j " r J Y j 

where, rj = ij' + rj\ 

Substituting Eq.(2.14) into Eq.(2.10) we obtain 

Qf = Qj mjCD 2 ( Y[ - ij' Y’[) 

Substituting Eq.(2. 17) into Eq.(2. 11) we can derive 

Mj R = + ij Q[ - mj co 2 rf Y[ - ( Ij - mjijrj 1 )co 2 Y' J L 

Collecting Eqs.(2. 16), (2.15), (2.17) and (2.18), and writing in 
form we obtain 


( Y l 

l R 

f Y l 

Y* 

l Q l 

II 

13 

2 

lY' 

Q 

l M ) 

! i 1 

M / 


where, the point matrix of the rigid body element, 

1 -rj 0 


[PM] j = 


0 


- mjGr 


1 


mi 


- m 


ljC0 2 Ij* 


ijCiAj' 
(Ij-mjrfijOco 2 


0 

1 


0 

0 

0 

1 


Substituting Eq.(2.7) into Eq.(2.19) we obtain 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

the matrix 
(2.19) 
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( 2 . 20 ) 


Y 

Y* 

Q 

M 




( Y\ r 

M 

l Q 

Mh 


where, the transfer matrix of the jth section, 



= [PM] j [FM] j = 


Oil Oi2 Ol3 Ol4 

O21 O22 O23 o# 

<X >31 O32 O33 O34 

_ O4I O42 O43 O44 


J jj-1 


The elements of the transfer matrix are listed below: 

<D n = ±<cospL + coshpL) + i-pr(sinpL - sinhpL) 

<J> 12 = —(sinpL + sinhPL) - jkcosPL + coshpL) 

2p 2 

0 13 = — 1— (sinpL - sinhpL) - (cospL - coshpL) 

2kp 3 2k P 

cD 14 = — 1 — (cosPL - coshpL) + (sinPL + sinhpL) 

2kp 2 2k P 


O21 = - ^(sinpL - sinhpL) 

<J>22 = ^(cospL + coshpL) 

023 = — - — (cosPL - coshpL) 

2kp 2 

024 = — —(sinpL + sinhpL) 

2k p 

<1)3! = - lmco 2 (cospL + coshpL) - ^Pmco 2 r L (sinpL - sinhpL) + ^P (sinpL + sinhpL) 

0 32 _ . -L m o) 2 (sinpL + sinhpL) + ^-mto^KcosPL + coshpL) - i-kp 2 (cosPL - coshpL) 
2p 2 

<D 33 _ . nK>ii( s i n pL - sinhpL) + m£^(cospL - coshpL) - kcosPL + coshpL) 

2kp 3 2kp 2 

<D 34 = . m(Q 2 ( cos pL - coshpL) - mi ^(sinpL + sinhpL) + |<sinpL - sinhpL) 

2kp 2 2k P 

<j) 41 = - — mco 2 r R (cosPL + coshpL) + y(I-mr R r L )o) 2 (sinPL - sinhpL) 

+ —kp 3 r( sinpL + sinhpL) - i-kp 2 (cospL - coshpL) 

2 ^ 
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$42 = - — moA^sinpL + sinhpL) - l(I-mr R r L )G) 2 (cospL + coshpL) 

2(3 ** 

- ^kp 2 r(cospL - coshpL) - lkP(sinpL - sinhpL) 

$43 = - — mco 2 r R (sinpL - sinhpL) - — L_ (I-mr R r L )o) 2 (cospL - coshpL) 

2kp 3 2kp 2 

- ^-(cospL + coshpL) - — (sinPL + sinhpL) 

1 2p 

$44 = - mw 2 r R (cospL - coshpL) + -J— (I-mr R r L )(o 2 (sinpL + sinhpL) 

2kp 2 2kp 

+ ^-pr(sinpL - sinhpL) - ^-(cospL + coshpL) 

3. Characteristic Equation; Eigenvalue and Eigenfunction 

After establishing the equation of motion of a global system by 
combining the transfer functions of all necessary elements, the natural 
frequency and mode shape function can be solved by satisfying the 
appropriate boundary conditions. 



Fig. 3.1 The Mathematical Model for NASA Mini -Mast 


As a mathematical model of the NASA Mini-Mast truss, it consists of two 
successive beam elements with two rigid bodies at the BaylO and Bayl8 
(Fig.3.1). Using the transfer function in Eq.(2.20) we can see 

G-p. 
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(3.1) 



have 


Now let us consider the boundary conditions. 

Y 0 = Y'o = 0 


At the fixed end we 
(3.2) 


at the free end, 


Q 2 = M 2 = 0 


(3.3) 


Applying the BC.'s to the Eq.(3.1) we get 


|v Y | 

R 

r A i 

0 1 
0 \ 

0 

= 2.0 < 

i Q 

\ 0 1 

2 

\ m) 


Rearranging the state vector we will have 



(3.5) 


The condition for Eq.(3.5) having non-trivial solution is that the determinant 
of the coefficient matrix equals to zero, that is. 


Det[A] = Det 


- 1 
0 
0 

_ 0 


0 <1>13 ^>14 

- 1 <t>23 *$24 

0 <t>33 <X>34 

0 C>43 <t>4i - 


(3.6) 


where, On's (i=l to 4, j=3,4) are the elements of the transfer matrix [<J>] 2 ,0- 
Eq.(3.6) is the so-called characteristic equation. Expanding the determinant 
in Eq.(3.6) we can rewrite the characteristic equation as 

O 33 O 44 - O 34 O 43 = 0 (3.7) 

Solving for the roots of the characteristic equation, Eq.(3.7), we can get the 
eigenvalue's (3's. To verify the theoretical derivation we have deduced the 
characteristic equations for two simple examples from the foregoing 
characteristic equation (Eq.3.7) as follows: 

For a cantilevered beam: 
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cosPLcosh(3L = -1 

For a cantilevered beam with a tip body; 


1 + 


cospLcoshpL = - pL ( cospLsinhpL 

Wb 


sinpixoshpL 


where W b is the weight of the beam; W is the weight of the body. These 
results are identical with those given in commonly used textbooks [8]. 

To establish the mode shape functions we must solve the Eqs.(3.5) and 
(2.19) simultaneously, that is, 


and 


- 1 0 

0 - 1 
0 0 

- 0 0 


$13 

$23 

033 

^43 



Y 

Y'l 

Q 

Ml 


= [PM] 




(3.8) 


(3.9) 


The later equation represents the compatible conditions of the deflections, 
forces, and moments between the two sides of the Bodyl. Because of the 
translational and rotary inertias of the Bodyl, the shear and bending 
moment have jumps at the connection point of the two beam elements, 
while the deflection functions are still continuous. 

According to the solution function, Eq.(2.4), we can express the state 
vectors in Eqs.(3.8) and (3.9) in terms of the coefficients A;, B<, C, and Di 
(j= 1,2) by using the following relations: J 


Q^ = A 1 k 1 p?-C 1 k 1 p? 
M^ = B 1 k 1 p?-D 1 k 1 p5 


Y V - AisinPiLi + BicospiL] + CisinhpiLi + DicoshPiLi 

Y j = AiPjcosPiLj - BiPisinPiLj + CiPjCoshpjLi + DiPisinhPiLi 
Qi" = ' Aikjp^cospiLj + BikjPisinpiLi + Qk^coshpiLi + DjkipJsinhpjLj 
= ‘ AjkipisinPiLi - BjkiPicosPiLi + CikjPisinhPiLi + DjkiPicoshPiLi 

Y ? = B2 + D2 
Y'f = A2P2 + C2P2 


(3.10) 
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Q? = A 2 k 2 p2 - C 2 k2p2 
M? = B 2 k 2 p2 - D 2 k 2 p2 


s Y 2 = A2sinp2L2 + B 2 COSP 2 L 2 + C2sinhP2L 2 + D 2 coshP2L2 
Y'| = Y 2 = A 2 p2COSp 2 L 2 - B 2 p2sinp 2 L 2 + C 2 p2Coshp 2 L 2 + D 2 P2sinhp 2 L 2 

Substituting the quantities in Eq.(3. 10) into Eqs.(3.8) and (3.9) we 
obtain a matrix equation in the coefficients Aj, Bj, Cj and Dj ( j= 1,2 ), 


kipi<J>i3 

klpfou - klPl<t>13 - klPl<l>14 

015 

Oi6 

a n 

klpl<t>23 

kiP?<t>» - kiPi<t>23 - kiP?«J>i* 

025 

026 

Ot? 

ki pl<I>33 

klp?<t>34 - klPl<t>33 - klP?<|)j4 

0 

0 

0 

klpl^>43 

klP?<t>44 - klPl<t>43 - klp?<l>44 

0 

0 

0 

Cts\ 

052 053 054 

0 

-1 

0 

061 

062 063 064 

-P2 

0 

-P2 

071 

072 O 73 O 74 

- k 2 p 2 

0 

k 2 p 2 

Ogi 

(X82 083 084 

0 

- k2p 2 

0 


CX 18 
a 2 s 
0 
0 
A 
0 
0 

k2p2 



A 2 

b 2 

C2 

'D2 


=[ 0 ] 


(3.11) 


The elements expressed in (Xij's in coefficient matrix are listed below. 


ais = - sinp 2 L 2 
ai6 = - cosp 2 L 2 
an = - sinhp 2 L 2 
= - coshp 2 L 2 


a 2 5 = - P2COSp 2 L 2 

a26 = p2sinp 2 L 2 
a 2 7 = - p 2 Coshp 2 L 2 
a 2 8 = - P2sinhp 2 L 2 


asi = sinPiLi - riPicosPiLi 
«52 = cosPiLi + ripisinpiLi 
a53 = sinhpiLi - riPicoshpiLi 
a 54 = coshpiLi - npisinhpiLi 

a6i = PicosPiLi 

a 6 2 = - PisinpiLi 

a 6 3 = PicoshpiLi 

a 64 = PisinhPiLi 
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CX71 = - miCO?sinpiLi + mito^p^ospiLi - kip^cosPiLi 
a 7 2 = - mi(o 2 cospiL] - miCo^pjsinPiLi + kjp^sinpiL! 
a 7 3 = - miO) 2 sinhpiLi + miW^PicoshpiLi + kip^coshp^i 
a 74 = - miG^coshpjLi + miW^pisinhpiLi + kipisinhPiLi 


a 8 i = - miC0 2 ifsinpiLi - (Ii-mu?rV)co 2 PicospiLi - rikiP^cospjLj - kip?sinpiLi 
a 82 = - mi0 2 ifcospiLi + (Ii-mjrf^^pjsinpjLi + nkip?sinp,Li - kip^cospiLi 
a 83 = - miCO^sinhpjLj - (Ii-mir?rV)co 2 p 1 coshp 1 L 1 + rjkjp^coshp^! + k^fsinhp^! 
a 84 = - mj^ifcoshPjLi - di-mjr^o^sinhp^ + nkip^sinhp^j + kiP^coshpjL, 

The solution to the Eq.(3.11) will give infinite sets of coefficients 
corresponding to each order of eigenvalue, which are usually called the 
modal participant coefficients. Assume that the solution is normalized with 
respect to D 2 , that is, D 2 =l. Then we will have a specific set of coefficients 

corresponding to an eigenvalue Pjj, which is now assumed in the form of 

Ai=ci i? Bi=C 2 i, € 1 = 03 ^ Dj = c 4i ForBeaml 

a 2 = C 5 i » B 2 = c 6i , C 2 = c 7i , D 2 =l For Beam2 (3 ' 12) 

Substituting the coefficients in Eq.(3. 12) into the solution equation Eq.(2 4) 
we obtain the eigenfunctions, or the mode shape functions, 

Y;(z) = / Cl i sin P 1 i z+C 2 . cos P]i z+c 3 i sinh Pi 1 z+C 4 ,coshp li z For Beaml 

\ c 5 i sinp 2 z+c 6 i cosp 2 z+c 7 sinhp 2 z+coshp 2 z For Beam2 ^ ^ 


4. Theoretical Response 

When the proportional damping is taken into account the Bernoulli- 
Euler beam equation will be 


d 2 y dy d\ 
m— - + c — + k — = 0 


dt 2 dt 


dz 4 


(4.1) 


where, c is a damping constant of proportionality which is assumed as 
c=2bm. After separation of variables and introducing the generalized 
coordinates Ti(t), we obtain 



[mYjYj] dz + 



[cYjYj] dz + 



[kY;Yj] dz = 0 
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which yields all zeros except for the one term in each when i-j according to 
the orthogonal property of the eigenfunctions, that is, 


mi Ti(t) + CiTi(t) + kiTi(t) = 0 


(4.2) 


where, 


mi = j mY|dz 

generalized mass 

J 0 

kj = k Yj dz 

generalized stiffness 

J 0 

Ci = c Y? dz = 2 bmi 

Jo 

generalized damping 


Expressing Eq.(4.2) in modal form we obtain 

Ti + 2^iC0i Ti +0)?Ti =0 (43) 

where. ^iWi-c^mi and o>i 2 = kM . The solution to Eq.(4.3) is now 

Ti ( t ) = e ( Aj cos 03<t t + Bj sin < 04 1 ) (4.4) 

where Ai and Bi are the coefficients dependent on the initial conditions. By 
superposition, the solution to the Eq.(4.1) can be written as 

(4.5) 


y(z, t) = X Y i( z ) e ■&* ( Aj cos w* t + B; sin co* t ) 


where, the eigenfunction Yi(z) has been derived in Eq.(3.13). Recall that we 
have defined that 


a 2 =k/m, 2b=c/m, and p 4 =co 2 /a 2 


(4.6) 


Then the damping ratio and the damped natural frequency co di can be 
expressed in terms of the parameters a and b, 


p. = c i — = _b_ and t0(j. = (0j V 1 - - V (api ) - b 2 

* ' 2miC0i ap 2 


(4.7) 


By superposition, finally, the solution to the Eq.(4.1) can also be written in 
terms of the parameters a and b, 
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y(z, t) = £ Yj(z) e bt ( A;cos t V (apf) 2 - b 2 + Bjsin t V (ap?) 2 * b 2 ) 


(4.8) 


5. Estimation of Modal Characteristics of NASA Mini-Mast 


Xhe maximum likelihood estimation method (MEE) is used to estimate 
the modal characteristics of the NASA Mini-Mast truss. The iterative formula 
for the MLE estimator has been derived in Ref.[9], 


9 = 0 O + 


I(V eyj ) T R-i(Ve yj ) 

Lj=i 


X ( V eyj ) T R- 1 ( yj - yjo ) 

Lj=i 


where, 

y ° nominal response calculated by using 0 o 
0 o nominal 0 vector 

^ gradient of y with respect to 0 
R covariance of the measurement noise 


(5.1) 


The unknown parameter vector 0 will be defined as 


0 = [a, b, A 1( A 2 , . . . , A n , B,, B 2 , . . . , BJ T 


(5.2) 


Because we have got the closed-form solution y(z,t) (Eq.4.8), the gradient 

V eyj can be easily obtained by directly taking the derivatives of y with 
respect to the unknowns. The closed-form expressions of the gradients 
have been derived, 


f< z ’t) - X Y i( z ) e ' bt -7==== [ - Ai sin t V (ap?) 2 -b 2 + B; cos t V (a|3?) 2 - b 2 ] 
da 1 V (ap?) - b 2 

£(*)- X Y i( z ) t e _ b 1 1 - Ai cos tV (ap?) 2 - b 2 - Bj sin tV (apf) 2 - b 2 

C/D j 

+ — b - [ Aj sin tV (apf) 2 - b 2 - B; cos tV (ap?)* - b 2 ] I 

V (aft?) 2 - b 2 1 

~ (z,t) = Y;(z) e - b [ cos tV (a(3f) 2 - b 2 
3Ai 

— — (z,t) = Yi(z) e - b 1 sin tV (apf) 2 - b 2 
3Bi 

(5.3) 
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Now we use Eq.(5.1) iteratively, considering yj as the measurements 

on a certain location of the beam at each time instant, and 7) 0 as the iterative 
response values calculated by using the updated 0 O , at the same location an 
instant When the innovation of the unknown parameter vector reaches the 

required criterion we may obtain the estimate 0 of 0. All the modal 
nrooerties are related to the parameters a and b (Eq.4.7), thus we can obtain 
thermo dal properties as long as these unknown parameters are determine . 

The test data is contained in Ref.[10], which was measured by one 
displacement sensor installed at Bay 18, mounted parallel to the flat face on 
the P corner joints of the structure and positioned to measure deflections 

norraaUo the face. estima ted frequencies which are compared with 

those obtained from Finite Element Analysis (FEA) and an Eigensystem 
Realization Algorithm (ERA) [11]. Fig.5.1 shows that the reconstructed 

response obtained from the estimated parameters and the measured response 
have a reasonably good fit 


6. Concluding Remarks 

This paper has demonstrated the principles for applying a t ran sfer 
matrix method to the parameter estimation of large space structures. The 
transfer matrix for thesystem with flexible beam elements and rigid bodies 
has been derived. The procedure for establishing natural frequency and 
mode shape has been described in detail. Maximum likelihood estimation 
method has served to conduct the parameter estimation. Comparing with 
me finite element model, the decrease in the number of unknown 
parameters by the present method is significant. The calculation, therefore, 
becomes highly efficient. The estimated results are compatible with those 

obtained by other traditional methods. ... w 

Further research is needed to formulate a more general method for 

more complicated structures. Some problems require 
transformation for non-perpendicular attachment elements. Transfer 
matrix for a branched structure must be considered. It is also desirable to 
develop a more efficient computer software based on the transfer matrix 
method, such as the new version of PDEMOD [12]. 
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Response (ft.) 


Table 5.1 Comparison of Estimated Bending Frequencies (tlx.) 
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Fig. 5.1 Comparison of Reconstructed and Measured Responses 
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